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ABSTRACT 

We examine the structure of adiabatically contracted Einasto profiles, using 
the prescriptions of Blumenthal et al. (1986) and Gnedin et al. (2004), and 
its impact on rotation curves. Adiabatically contracted halos display a central 
power index of ~ 0.7 ±0.1 for nearly all values of the Einasto shape parameter 
a, and are well fit inside < 1 kpc by double power laws. However, attempts to 
fit exponential disc and uncontracted halos to adiabatically contracted rotation 
curves yield disc masses and central power indices that are too large. We also 
determine whether or not the rotation curve of NGC 6503 displays evidence of 
adiabatic contraction using previously published bar formation constraints, and 
find that NGC 6503 most likely has a minimally contracted halo. However, this 
conclusion depends on the correct choice of circular velocity curve. 

Subject headings: galaxies: individual (NGC 6503) — galaxies: kinematics and 
dynamics — galaxies: halos 



1. Introduction 

One of the most important unsolved problems of cosmology and galaxy evolution is how 
baryons affect dark matter halos. Over the past two decades, high resolution cosmological 
simulations have effectively settled the question of what form dark matter halos take in the 
ACDM cosmogony, which has been very successful in accounting for the large scale structure 
of the universe. However, on galactic scales, detailed studies of the structure of galaxy halos, 
as inferred from rotation curves and lensing measurements, have produced mixed results, and 
it is critical to understand how baryons affect dark matter halos to resolve the discrepancy. 

Traditionally, d o uble power laws have been used to fit cosmological iV— body halos 



( jNavarro et al.l 119971 : iMoore et aJj Il999l : iDiemand et al.l 120051 ) , the most general form for 
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which is given by 



p(r) 



(r/r s ) 7 1 + (r/r B y 



where p s and r s are the scale density and radius respectively, and the triplet (/3, 7, 5) governs 
the shape of the log-log profile. Commonly, (3 and 5 are set to 1 and 3 respectively, allowing 
only the central power index to vary: 



p{r) 



Ps 



(r/r s ) 7 (1 + r/r s ) 



3-7" 



(2) 



We call Eq. [2] the N FW-type profile. Setting 7 = 1 produces the traditional NFW formula 
(INavarro et al.lll997l ). and the resulting two parameter profiles were thought to be universal. 
Recently, high resolution simulations have sugge sted that simulated halos vary slightly i n 
shape, thus requiring an extra sh ape parameter (IMerritt et al.ll2006l : INavarro et al.ll2010l ). 
The adopted profile is usually the lEinastol (119651 ) profile, given by 
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This form is identical to that used by INavarro et al.l ( 120101 ) ; and ph are the radius and 
density of the peak of the r 2 p profile, respectively, and a is the shape parameter. Eq. [3] is 
formally identical to the classical Sersic (1968) 1/n profile that is often used to model the 
projected surface brightnesses of early type galaxies and bulges; however, the Sersic profile 
is applied to the projected radius, while the Einasto profile is applied to the spatial radius. 
For studies of halo structure, Eq. [3] is preferable because the constants are directly related 
to the easily calculated r 2 p quantity. 

Modifications to Eq. [3] will occur because of baryon condensation as galaxies form. The 
simplest prescription for baryon-affected dark matter halos is adiabatic contraction, which 
posits that halo orbits are purely spherical and do not cross, and that ba ryon condensation 
occur s adiabatically, and therefore that angular momentum is conserved (IBlumenthal et al. 
19861 ). The adiabatic contraction condition may be expressed as 



n M{n) = n [M f (r f ) + M b (r f )] 



(4) 



where M(r{) is the initial mass of dark matter and baryons at radius n, Mp(rf) is the final 
mass of dark matter at radius rf, and Mb(rf) is the mass of the final baryon configuration 
at rf (usually Mb is an exponential disc, but may in principle be any distribution). The 
equation may be easily solved for r f using iterative methods. Adiabatic contraction leads to 
a pronounced concentration of mass in the centre of a halo, steepening the cusp. 
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While simple, adiabatic contraction is unlikely to provide a complete description of 
how baryons impact dark matter halos, because it does not account for radial motions 
(essentially assuming that the radial action is zer o). Additionally, s i mulations suggest 



that angular mo mentum is not, in fact, conserved (IGnedin et al.l 12004 ; iDuffy et al.l 12010 



Book et al. 
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Gnedin et al- 



ii). M odified adiabatic contraction prescription s, such as those introduced by 



(12004). a nd simulations, such as those done by ISeUwood fc McGaughl (120051 ) 



and lAbadi et al.l (l2010f l. still predict a relative contraction of matter in the centre, but not 
as much as Blumenthal's original pre scription because the radial motions function as a pres- 
sure support resisting compression. IGnedin et al.l ( 120041 ) modify the adiabatic contraction 



prescription by replacing r in Eq. HJwith r, the orbit-averaged radius, as follows: 

nM (n) = r t [M f (f f ) + M b (f f )] 



(5) 



The relationship between r and f is given approximately by x = Ax w , where x = r/r v - ir , 
A = 0.8 and w = 0.8. This substitution i s desi gned to account for the fact that halo orbits 
are generally not spherical. lAbadi et al.l (120101 ) find a relationship given by 
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where B = 0.3, c = 2 and Mt/M x = Mfc) /\M r (r f ) + MJ r f )}\ this is, however, a statisti- 
cal average of severa l simulations. The Abadi et al.l (120101 ) prescription produces still less 



contraction than the IGnedin et al.l (120041 ) prescription. 



The full effect of baryons on dark matter halos remains an open question, as multi- 
ple effects must be accounted for to obtain a full picture of this process. There is am- 
ple evidence from A^— body simulations that t he contraction can probably be reversed by 



supernova feedback in the centres of galaxies (jNavarro et al. 1996 



Read fc Gilmore 2005 



Mashchenko et alJbood : iGovernato et al.lboioh . while other auth ors propose that the con 



traction may be revers ed or reduced by infalling gaseous clumps (lEl-Zant et al 



Elmegreen et al.ll2008f). by the effect of a bar torquing the halo (IWeinberg fc Katz 



2001 



2004 



2002 



Hollev-Bockelmann et al.l 120051 ) . by preprocessing of primordial halo clum ps (|Mo fc Mao 
20041 ) . or by the combined effect of dynamical friction and angular momentum (iDel Popolo fc Kroupal 



20091 ) . These factors complicate attempts to use adiabatic contraction and related theories to 
model observed rotation curves. Indeed, most measurements of galaxy halo profiles conclude 
that LSB galaxies have much more shallow cusps than implied by cosmological simulation s 
jFlores fc Primacklll994l : Ide Blok et aDboOll ; Isimon et allbood : IrCuzio de Narav et alil20061 ) , 
suggesting that adiabatic contraction is reversed; note that adiabatic contraction always pro- 
duces cuspy halos, even if the unc ontracted haloes are cored (cf. Fig. 2 of Blumenthal et al. 
1986 and Dutton 2005). However, lOh et al.l (1201 ll ) find that they can reproduce the obser- 
vations of LSB galaxies by including star formation, supernova feedback and extragalactic 



3 



UV heating, all of which overwhelm the effect of adiabatic contraction. Meanwhile, several 
studies of HSB galaxies suggest that th e data is consistent with cuspy and cored profiles, even 



with multiple constraint s factored in (IDutton et al.l l2005t Ide Blok et al.l 120081 ; lAuger et al. 
2f)10l : ISchulz et alibmnh . 



There is moreover no definitive proof on whether adiabatic contra ction is a val i d gen- 
eral description of baryon affected dark matter halos. For disc galaxies, Kassin et al.l (120061 ) 
find that adiabati c contraction yields poorer fits to the rotation curve, especially if the disc 
dominates, while IDutton et al.l (120051 ) find rotation curves with s ubma ximal discs equally 
well fit by both contracted and uncontracted halos. ISchulz et ah] (120101 ) find that adiabat- 
ically contracted NFW profiles provide an excellent fit for their sample elliptical galaxies, 
while lAuger et al.l (120101 ) find that adiabatic contraction significantly overestimates the dark 
matter content of early type galaxies. The discrepant conclusions reflect a broader degener- 



acy b etween halo contraction and the choice of IMF (jNapolitano et al.l l2010t IDutton et al. 



201ll ). but the very real possibility exists that the effects of baryons change significantly be- 
tween different classes of ga l axies, and even between d ifferen t galaxies of the same class: the 
simulations of lAbadi et al.l ( 120101 ) and iTissera et al.l (120101 ) suggest that the halo assembly 
history, which is certainly not identical for all galaxies, plays a crucial role in determining 
the halo response to the baryons. 

In this paper, we restrict our analysis to how Einasto profiles and their associated 
rotation curves are affected by the adiabatic contraction prescriptions of Blumenthal et al. 
and Gnedin et al. In particular, we would like to know if it is possible to unambiguously 
identify adiabatically contracted rotation curves in disc galaxies, or if sufficient uncertainty 
exists that this is not possible. The structure of this paper is as follows. In Section 2, we 
discuss how density profiles and rotation curves are affected by adiabatic contraction, using a 
model halo to investigate these issues. In Section 3, we use published rotation curve data for 
NGC 6503 to determine if the dark matter content can be constrained. Section 4 discusses 
and summarizes our conclusions. 



2. The Structure of Adiabatically Contracted Halos and Rotation Curves 

There are three parameters characterizing adiabatically contracted (AC) halos: the 
baryon fraction f^, the baryon scale radius (which is the exponential scale length in our 
case), and the radius out to which baryons condense, r outer . This last radius need not in 
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principle be the virial radius r vir £] To examine the structure of AC halos, we begin with 
a model Einasto profile (Eq. [3]) with = 22.65, ph = 0.00058265 and a = - 17. T hese 



parameters are chosen to match the mass-concentration relation of iMaccio et al.l (120081 ) and 
to produce a halo and disc combination whose rotation curve peaks at ~220 km s _1 . For 
comparison, we include a halo with identical ph and with a = 0.5. We then apply the 
Blumenthal and Gnedin adiabatic contraction prescriptions at varying degrees of contraction 
given by f£ = 0.01, f£ = 0.075, and f£ = 0.15. We use GalactlCS units for this exercise 



([Widrow et al.ll2008l ). in which G = 1, the unit of length is 1 kpc and the unit of velocity is 



100 km s 1 , so that the unit of mass is equal to 2.325 x 10 9 M Q . 

Fig. [1] shows how the r 2 p profile of these AC halos compares to standard Einasto profiles. 
At higher degrees of contraction, AC profiles display a double inflection just beyond the peak 
of the r 2 p profile that neither Eq. [TJ nor Eq. [3] can reproduce. The double inflection is not 
apparent on a logp — logr plot, but is very clear on a log r 2 p plot. Adiabatically contracted 
profiles also display a nearly linear structure inside the r 2 p peak, indicating power law 
behaviour. Einasto profiles cannot reproduce the power law behaviour, but NFW-type and 
double power laws exhibit power law behaviour as r — > 0. 

To determine if NFW-type or double power laws can fit the inner part of an AC halo, we 
plot fits to AC halos in Fig. [2] using Eq. [2] and in Fig [3] using Eq. [TJ For both prescriptions, 
the NFW-type profile fails to provide an effective fit, but the double power law yields a much 
better fit, provided ~ 0.5 — 1 for a = 0.17 and /3 ~ 1.5 — 2 for a — 0.5. The fits become 
inferior as ft, decreases (that is, in cases where the final halo is close to the original) at larger 
a, but are otherwise very good. Thus, the central power 7 C obtained by direct calculation is 
identical to 7 in Eq. [TJ 

The values taken by 7 display a remarkable consistency. Fig. [TJ provides a scatter plot 
showing how the central power 7 C varies with a, f^, r b and r b /r h . At a > 0.3, 7 ~ 0.8 
as r — 0. For a < 0.3, 7 displays larger variations and becomes larger as a — > 0. The 
consistency of 7 is maintained across the other three parameters in the figure, including 
cases for which 0.1 < a < 0.2, suggesting that 7 is not strongly dependent on any of these 
parameters. 

The consistency of 7 suggests a possible test for adiabatic contraction. Double power 
law fits to the rotation curve inside the peak of the r 2 p profile could be used to determine 7 
and thus if a given disc galaxy displays evidence of adiabatic contraction. The average 7 for 
the Blumenthal and Gnedin prescriptions is 0.75 ± 0.17 and 0.66 ± 0.08, respectively. For 



1 We define the virial radius as the radius within which the mean density of the halo is 200 times the 
critical density p cr i t . 
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cases where 0.1 < a < 0.2, the Blumenthal and Gnedin prescriptions yield 7 = 0.75 ± 0.17 
and 0.61 ± 0.08 respectively. Although there is considerable overlap in the two values for 
7, Blumenthal's prescription yields larger 7 with more scatter, suggesting that it may be 
possible to distinguish between the two types of adiabatic contraction prescriptions as well. 
Using Eq. [1] rather than Eq. [2] is critical, however, because on observationally accessible 
scales 7 is not constant, and because Eq. [2] can only provide an effective 7 over the range of 
radii probed, which limits its usefulness as a diagnostic parameter. Eq. [T] yields the true 7 
because the variation in power over a given range of radii can instead be captured by the /3 
parameter. 

However, the main probe of halo structure is rotation curves, which depend on the total 
mass (i.e., the integral of p) and hence the double inflection, the subtle differences between 
NFW-type and double power laws, and the differences between the Blumenthal and Gnedin 
prescriptions may be difficult to detect in real galaxies. Fig. [5] displays model rotation curves 
for the Einasto profile found in the left panel of Fig. [TJ adiabatically contracted at different 
values of f^. To to further examine this issue, we now attempt to fit pure Einasto and 
NFW-type profiles to a mock rotation curve for an AC halo, and apply it to a real galaxy 
in the section following. 



2.1. Rotation Curve Fits to AC Halos 

There have bee n relatively few attem pts in the literature to directly fit Einasto profiles 



to rotation curves - 



Chemin et al.l (120091 ) have done so for M31, but they were not able to 
constrain the halo profile. Moreover, many attempts to fit any kind of density profile to 
rotation curves ignore the effects of baryon condensation. In the case of LSB galaxies, this 
approach is probably justified because of the small disc mass. There have been far more 
attempts to directly fit NFW-type profiles to rotation curves. It is however interesting to 
determine how well AC halos can be fit by pure Einasto profiles and double power laws. 

There are three parameters governing the Einasto profile, and three governing adiabatic 
contraction as indicated in §2.1. These six parameters are summarized in Table [H It should 
be noted that r ou t er is not determined by the observed extent of the rotation curve; r outcr 
impacts the total mass of baryons, and therefore differences in this parameter will impact 
the entire rotation curve, even if all other parameters are kept constant. 

We adopt the same Einasto halo as in §2.1 for consistency with the concentration- 
mass relation, and include a disc of a fixed mass mdi SC = 12 and scale length = 3. The 
resulting rotation curve peaks at about 220 km s -1 . We then apply the adiabatic contraction 



(3 



prescription, using the same halo parameters, at different values of /b : fb — 0.025, 0.05, 
0.075 and 0.1. We reproduce the initial disc mass by adjusting r outer , so that the disc mass is 
identical for all /b- The resulting model rotation curves are sampled linearly in radius with 



error bars and resolution roughly consistent with THINGS data published in Ide Blok et al. 



( 120081 ) . to which noise is added. We then attempt to fit these rotation curves using pure 
Einasto profiles, holding the disc scale length fixed since this quantity may be constrained by 
surface brightness profiles. We use a Bayesian/MCMC approach to find the best fit for each 
case, and then carry out the same procedure for Gnedin's adiabatic contraction prescription. 
The original and AC rotation curves are shown in Fig. [5] for reference. 

The resulting fits are found in Fig. [6] and the posterior mean and error bars are found 
in Table El with GalactlCS masses converted back to physical masses. The fits are excellent, 
suggesting that AC halos are well fit by Einasto profiles. However, the fits produce smaller 
a than found in the original halo, and the disc mass required to achieve these fits is about 
40% larger than used to construct the rotation curves in each case, suggesting that attempts 
to fit AC rotation curves with Einasto profiles overestimate the disc mass. This is because 
AC rotation curves have steeper slopes near the centre, and heavier discs are required to 
reproduce this feature when using a pure Einasto halo. Using Gnedin's adiabatic contraction 
prescription alleviates the problem somewhat, yielding marginally lighter discs. 

While the best fits overestimate the disc mass, it may be possible to obtain good fits 
by fixing the disc mass to its correct value|§ In this case, we find that the fits are not as 
good but improve as /b increases and improve when Gnedin's prescription is used as well; 
in general, however, these fits are inferior to those that allow the disc mass to vary and are 
not consistent across various values of the adiabatic contraction parameters. 



2.2. Fits Employing NFW-type and Double Power Profiles 

The fact that AC Einasto profiles exhibit power law behaviour inside suggests that 
Eqs. [T]and[2]may produce better fits to the rotation curves of AC halos. We first carry out fits 
with Eq. [2] in which (3 and 5 are fixed to 1 and 3, while the central power index 7 is allowed to 
vary. These fits are otherwise identical to the fits carried out in the previous subsection. The 
resulting fits are found in Fig. [7] and are essentially identical to the Einasto fits. However, 
they too overestimate the disc mass, in some cases by more than 50%. Allowing /3 and S 
to vary still overestimates the disc mass (Fig. [8]). We conclude that NFW-type and double 



2 Note that the error bars in Table[2]are formal 1—a error bars obtained from MCMC. Qualitatively, such 
error bars may not be large enough to provide an acceptable range of values. 
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power profiles offer similarly good fits to the rotation curves of AC Einasto halos, although 
disc masses are overestimated. 

Note that 7 is found to be ~1.2, which exceeds the value of ~ 0.6 — 0.8 for these halos, 
as indicated in Figs. [2] and [3j We therefore also attempt to fit a double power profile with 7 
fixed to 0.8, finding that the fits require very massive discs (more than twice the true mass). 
Moreover, attempts to fit double power laws to these rotation curves with a fixed disc of the 
correct mass fail to yield the correct value of 7. We conclude that although AC profiles may 
be described by a power law inside a certain radius, NFW-type and double power profiles 
cannot be used to fit the full rotation curves of AC halos because these fits overestimate the 
disc mass, and because the fitted 7 exceeds the true 7. 

2.3. Fitting only the Inner Rotation Curve 

While using the full rotation curve yields incorrect values for 7, fits that only use the 
inner rotation curve may produce better values. Moreover, we would like to determine if it 
is possible to distinguish between the Blumenthal and Gnedin prescriptions using rotation 
curve data. Therefore, we carried out the same fitting procedure as in the previous section 
for double power fits, this time truncating the rotation curve at 5 kpc (well inside the peak 
of the r 2 p profile for all halos). 

The resulting values for /3,7, and S are found in Table [3j These may be compared with 
the 'true' values found in Table HJ which were obtained by fitting the r 2 p profile directly. 
We find that 7 < 0.9 for all fits; in most cases, the true 7 values are not reproduced by 
the fits, although these estimates are much better than those found in the previous section. 
However, the error bars are quite large (especially for j3) and the correct values are generally 
captured at the 2 — a level. In addition, there is considerable overlap between the fitted 
values for the Blumenthal and Gnedin prescriptions. Thus, not surprisingly, the mock data 
cannot distinguish between the two prescriptions, probably because the error bars are not 
small enough. We conclude that it is essentially impossible to distinguish between the 
Blumenthal and Gnedin adiabatic contraction prescriptions unless the observational error 
bars, and possibly resolution, improve substantially. 

3. NGC 6503 as a Test Candidate 

NGC 6503 is an attractive candidate to test how baryons affect dark matter halos. It 
is mostly isolated, so the impact of mergers and tidal interactions is limited, and therefore a 
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direct comparison of the simple adiabatic contraction prescription is more ap plicable to this 



galaxy than other galaxies. Published model circular velocity curves exist from lPugrielh et al. 



(120101 ) (hereafter PWC), which will allow a more accurate comparison of adiabatic contrac- 
tion with the true mass distribution of the halo. Moreover, the fact that NGC 6503 is clearly 
strongly cusped invites a direct test of the adiabatic contraction hypothesis. 

The appropriate choice of rotation data must be made. Ideally, we would use data that 
perfectly traces the gravitational potential, and subtract the effect of the disc and bulge. 
However, rotation curve measurements consist of ionized gas, H I, or stellar kinematics, 
none of which trace the potential as precisely as we require. Stellar rotation suffers from 
asymmetric dri ft (cf. Eq. 4.228 of B inney & Tremaine 2008), while H I has dispersions of 



~5-25 km s _1 (ITamburro et al.ll2009i ). As we argue in PWC, the best method of obtaining 
the circular velocity may be to calculate the asymmetric drift of the stars and add that to 
the observed stellar rotation. Thus, we will use the circular velocity from PWC as part of our 
data set. Since the asymmetric drift is itself highly sensitive to several effects (in particular, 
it can be made very close to zero if one assumes that the velocity ellipsoid is spherically 
aligned), we will also test adiabatic contraction assuming that the gas rotation does trace 
the gravitational potential. This scenario should be treated as a lower limit on the possible 
circular velocity curves. Fig. [9] presents both data sets. 



3.1. Potential complications 

The effect of a bar on the rotatio n curve is a potentia l complication, as NGC 6503 is 



now known to possess an end-on bar (IFreeland et al.ll2010l ). Because the bar is end-on, it 



is unlikely to significantly affect measurements of either the stellar or gas rotation, and we 
therefore do not anticipate any significant bias to our results. The presence of the bar is 
intimately connected to the question of the correct disc density profile to use. PWC found 
that an inner truncated disc may provide a good fit to the observati ons (their scenario K) . 



The close relationship between measured inner truncations and bars ([Anderson et al.ll2004l ) 
suggests that the inner truncation in NGC 6503 may be due to the bar, reducing the problem 
to that of fitting an exponential disc and evaluating the consistency of the resulting fits with 
the disc stability properties required to reproduce the observed bar. To evaluate disc stability, 
we determine the stability parameters Q and X, which measure local and global stability, 
respectively. 

A more serious complication is the possible presence of a massive bulge. The bulge 
of NGC 6503, like most other late-type galaxies, is a pseudobulge, and is therefore likely 
to have been formed by secular evolution. In such a process, bars and spirals cause gas 
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to lose angular momentum, funnelling the gas into the centre of the galaxy where stars 
form. Therefore, the conditions for adiabatic contraction are not satisfied and we restrict 
the analysis to cases where baryon condensation only results in an exponential disc. However, 
we must still account for the possibility of a massive bulge impacting the observed rotation 
curve - PWC have shown that the bulge mass is < 5% of the disc mass, arguably too large 
to be ignored. We account for this by adding the circular velocity due to the bulge found in 
PWC to the model rotation used in our fits. 



3.2. Our tests 



We first note that direct measurements of the density profile using the rotation curve 
are infeasible. In Fig. [10] we show the r 2 p profile for NGC 6503 obtained from the circular 
velocity and ionized gas data. The profile obtained from the circular velocity curve is shown 
assuming no disc and assuming a disc of 3.5 x 1O 9 M . The profile displays two peaks (owing 
to the kink in the circular velocity curve at ~ 7 kpc), but there is a clear global maximum at 
2.5 kpc (assuming no disc) and 10 kpc (assuming a massive disc). The central power cannot 
be well constrained using the available data, however. The r 2 p profile for the ionized gas 
data appears as little more than random noise. 

We therefore adopt the following approach. We employ two separate data sets for two 
different suites of models: the circular velocity fro m PWC fo r which the suite of fits is labelled 
C, and the H/3/H I data from lBottemal (119891 ) and iBegemanl (119871 ). for which the suite of fits 
is labelled I. The baryonic components may include a bulge and disc, or solely a disc. There 
are thus four possible ways to conduct the fit. Each set of four runs is impl emented assuming 



three contraction prescriptions: the adiabatic contraction prescription of iBlumenthal et al. 
(11986 ) which we label B, the modified adia batic contra c tion p rescription of iGnedin et al. 



(120041 ). labelled G, and the prescription of lAbadi et al. (boiok labelled A. Each of these 
prescriptions is applied to the standard Einasto profile. Finally, we conduct a final suite of 
fits assuming no baryon-induced changes to the standard Einasto profile, which we label E. 

In principle, we may fit either the total rotation curve or the halo rotation obtained by 
subtracting the baryonic rotation from the total rotation. We have done both for each fit 
described above and verified that the results are consistent with each other, so for clarity we 
only present the results from fitting the total rotation curve. 

We use a Bayesian/MCMC approach to find the best fit for each case. The choice of 
priors may b e either uni f orm, o r derived from the c once ntration- mass relati on obtained by, 
for example, iNeto et all (J2003), iMaccio et all j2008h . or iKlvpin et all ibuioh . Each of these 
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authors finds a slightly different relation, but all find a tightly constrained concentration- 
mass relation of the form C = kM^ h where k ~ 10 and e ~ —0.1, severely restri c ting the 



available range of parameters. In this work the C-M relation from iMaccio et al.l ( 120081 ) is 
used; because C and M vir depend on ph, r h, and a, the C-M relation implicitly defines 
priors for these parameters. We find little differences in our results for both sets of priors, 
so in what follows we present results for the C-M priors only. For the remaini ng; input 



param eters, we assume uniform priors. Although constraints on exist from iBottema 



( 119891 ) and PWC, we allow this parameter to freely vary to determine how well an adiabatic 



contraction prescription can reproduce these constraints. 

Two criteria must be used to constrain the adiabatic contraction process: first, iden- 
tifying which model produces the best fit; second, determining consistency with the bar 
formation constraints found in PWC. These constraints require that the disc not be vio- 
lently unstable, but because NGC 6503 is barred, a sufficiently large disc mass is required. 
PWC reported a best fit disc mass of 3 ± 0.4 x 1O 9 M for models whose circular velocity 
curve matches the one employed here. Only the upper end of this range is valid because lower 
mass discs do not form a bar. A better proxy for disc stability is the X parameter, because 
there is a strong correlation between the strength of the resulting bar and the X value of the 
disc. PWC found that (X) ~ 2.7 was required to produce a bar of the type found in NGC 
6503. Higher (X) values failed to produce a sufficiently strong bar; as a general cutoff, we 
assume that (X) < 3 is required to form a bar. Models with low (X) are likely to be too 
violently unstable to properly model the galaxy, but the cutoff is less well-defined. The Q 
value provides a better-defined lower bound on stability: It is critical that the Q values of 
the discs we find are sufficiently large to avoid fragmentation instability, which may occur if 
Q < 1, and so we also require that Q > 1 at all radii. 



3.3. Results 

Representative fits for each scenario are found in Figs.[TT]and[T21 along with the standard 
X 2 statistic. The x 2 mav be compared between models from suite I and between models from 
suite C, but should not be compared across the two types of fits because the error bars for 
suite C are not based on observational data. These error bars were obtained by randomly 
selecting several models from PWC and manually evaluating the standard deviation of the 
circular velocity values at each radius. The x 2 show relatively little variation, suggesting 
roughly equally good fits across all models; the lowest x 2 are obtained for fits to the Abadi 
prescription and to uncontracted halos, suggesting that the degree of contraction in NGC 
6503 may be small. Interestingly, the bulgeless fits are somewhat better than the bulge fits 
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for suite I, but there is little difference between the two for suite C. 

The factor producing the largest difference between models is the rotation curve used 
for the fits. In particular, fits from suite I invariably produce disc masses that are much 
larger than the fits from suite C. The baryon fraction for suite I is ~ 0.07 and the halo 
masses are about 7.5 x 1O 1O M , producing disc masses of ~ 6 — 7 x 1O 9 M . By contrast, fits 
from suite C yield /b ~ 0.03 and m-disc ~ 2.5 — 3.5 x 10 9 M Q . As Table [5] shows, disc masses 
tend to increase as the effect of baryon condensation decreases, i.e., pure AC halos produce 
the lowest mass discs and uncontracted halos produce the heaviest discs. 

To assess the bar stability of these discs, we can compare the disc masses for suite C 
directly to the results of PWC. The standard adiabatic contraction scenario produces discs 
that are too light and cannot form a bar, but the other scenarios produce discs that may 
be consistent with bar formation, although here too disc masses may be somewhat low. A 
similar comparison is not possible for suite I, so we instead plot Q and X on Fig. [13] for each 
fit, and present (X) in Table |5j The figure shows that, except for the standard adiabatic 
contraction scenario, Q is very close to or below 1. Below 1, the dominant instability is 
disc fragmentation and violent instability may result. On this basis, we reject most of the 
models from suite I. For the standard AC fits to the ionized gas data, the disc masses are 
systematically slightly smaller and hence Q and X are slightly larger, consistent with the Q 
and X values found in PWC that produced bars but are not violently unstable. Therefore, 
a scenario in which pure adiabatic contraction occurs and the circular velocity is accurately 
represented by the ionized gas curve is consistent with the data, but a scenario in which 
adiabatic contraction causes the circular velocity curve found in PWC is not consistent with 
the data because the discs of these models are too light (and hence X too high) to form 
bars. The reverse holds for scenarios G, A, and E, because the disc masses in those scenarios 
are higher and lead to lower Q and X. From PWC the bar formation cutoff occurs around 
Amm ~ 2.2 and therefore the circular velocity fits for scenarios A and E may be consistent 
with the observed bar. There is sufficient ambiguity in the derived masses and the derived 
X and Q values, however, that we are reluctant to make a definitive statement about the 
degree of dark matter contraction and about which models are most consistent with the 
data. 

Thus, the overall picture that emerges suggests that stronger halo contraction correlates 
with lighter discs. This is because the rotation curve is fixed, and so stronger contraction 
means that the halo contributes more to the inner rotation curve. This trend suggests 
two possibilities. First, the circular velocity curve found in PWC is correct and the halo 
contraction is not as strong as predicted by the standard adiabatic contraction prescription. 
Second, the correct circular velocity curve lies between the one found in PWC and the ionized 
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gas curve; in this case, disc masses are higher and the halo must be more strongly contracted 
relative to the first scenario. The two scenarios may be visualized by plotting a proxy for 
the degree of contraction and a proxy for disc stability. The stability proxy is simple: the 
best measure of disc stability is the X parameter. The contraction proxy is trickier. The 
natural choice would be the ratio of initial to final concentrations, but the concentration is 
defined in terms of the virial radius, and our approach assumes that the virial radius and the 
outer radius for baryons r oute r need not coincide. Therefore, we adopt an alternate approach 
in which the contraction proxy is given by the ratio of halo mass inside after and before 
contraction, R = M f (r b ) / Mi(r b ) . 

These two quantities are plotted in Fig. [HI demonstrating how higher contractions are 
correlated with higher (X), and also reveals that most of the models are located outside 
the ideal stability range to produce weak bars - in general, the fits to the circular velocity 
from PWC are too stable unless we assume almost no adiabatic contraction, while the fits to 
the ionized gas have such low (X) that they are likely to be too unstable unless we assume 
maximal adiabatic contraction. The ionized gas fits, represented as triangles on this figure, 
form a lower limit for X because the true circular velocity cannot lie below the ionized 
gas curve. It is possible to picture suites of models that employ differing circular velocity 
curves; on Fig. [TU these models would lie between models I and C, falling inside the ideal 
stability region. Therefore, the most likely scenario for NGC 6503 is that the circular velocity 
curve lies between the circular velocity determined in PWC and the ionized gas curve. If 
we instead assume that the circular velocity from PWC is correct, models that employ a 
standard Einasto profile without modification are most consistent with the desired stability 
properties. The shape parameter a is lower than found in the other models at a ~ 0.10, 
placing it somewhat below the values obtained from theory. 

It is interesting to note that the shape parameter a is also severely affected by the choice 
of rotation curve to fit, as seen in Fig.[T3 Suite C yields a ~ 0.1—0.15, n early consistent with 



(albe it slightly smaller than) those found in cosmological simulations ([Navarro et al.l 12004 



20101 ). For suite I, we find a ~ 0.6 — 1.0, wildly divergent from theory, and moreover that the 
PDFs are not well constrained. This striking difference testifies to the sensitive dependence 
of halo shape on the rotation curve used for the fits, and underscores the difficulty of inferring 
halo parameters from rotation curve observations. 

The choice of whether or not to include the bulge in the fits generally makes little 
difference to the overall result, because the bulge is small enough that its effect can essentially 
be replicated by altering the disc scale length slightly. As Fig. [TBI shows, bulgeless models for 
suite I produce discs that are more centrally concentrated. For suite C, the bulgeless models 
produce discs of slightly larger r b . Apart from this detail, however, the resulting fits are not 
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much different and the disc scale lengths produced by our fits are consistent with the value 
determined in PWC of 1.3 kpc, or nearly so to within ~10 % (Tabled]). 

4. Conclusion 

We have investigated the structure of adiabatically contracted Einasto profiles for the 
case where the baryons form an exponential disc, and obtained a number of interesting 
results. 

First, AC halos display power law behaviour as r — > and are well fit inside the peak 
of the r 2 p profile by double power laws, even in cases where the contraction is minimal. The 
central power index 7 is given by 0.75 ± 0.18 for Blumenthal's prescription and 0.61 ± 0.08 
for Gnedin's prescription for cosmologically motivated values of the shape parameter a. 
This result suggests a test for adiabatic contraction for various galaxy populations, on the 
assumption that uncontracted dark matter halos are well represented by Einasto profiles. 

Second, attempts to fit AC Einasto profiles and their exponential discs using uncon- 
tracted halos overestimate the disc mass significantly, suggesting that adiabatic contrac- 
tion must be accounted for when fitting mass models to rotation curves. However, the fits 
themselves are good for Einasto, NFW-type, and double power profiles. Double power and 
NFW-type fits using uncontracted halos tend to overestimate 7; better estimates for 7 are 
obtained when fitting only the inner rotation curve. 

Third, we have been able to fit the circular velocity curve of NGC 6503 to three different 
adiabatic contraction prescriptions, and found that obtaining the correct prescription is 
highly dependent on the correct choice of circular velocity curve. Assuming the rotation 
curve found in PWC, there is almost no contraction to the halo, but lower circular velocity 
curves may imply stronger contraction, up to the magnitude implied by Blumenthal et al.'s 
adiabatic contraction prescription. Thus, the available data fails to properly constrain the 
degree of halo contraction, underscoring the need for accurate determination of the circular 
velocity curve of spiral galaxies. 

The author wishes to thank L. M. Widrow for helpful discussions and critical com- 
ments on earlier versions of the manuscript, and the anonymous referee for comments and 
recommendations that greatly improved the manuscript. 
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Fig. 1. — Examples of adiabatic contraction applied to two different Einasto profiles, at 
differing degrees of contraction. The top panels use Blumenthal's prescription, while the 
bottom panels use Gnedin's prescription. For both halos, ph = 0.00058265, r h = 22.65; for 
the halo on the left, a = 0.17, while the right panels use a halo for which a = 0.5. The 
contracted halos are given by rb — 3, f£ — 0.01, by rb = 3, f£ — 0.075, and by rb = 1, 
/b = 0.15, in order from smallest to largest degree of contraction, for all four panels. Green 
identifies the original halos and black identifies the contracted halos. The contracted halos 
display a relatively constant inner power and, in the case of stronger contraction, a double 
inflection near or just outside the peak. Einasto profiles cannot reproduce either of these 
features. 
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Fig. 2. — Fits to two different Einasto profiles at differing degrees of contraction, using NFW- 
type profiles (Eq. |2J). The top panels use Blumenthal's prescription, while the bottom panels 
use Gnedin's prescription. The profiles and colours are as found in Fig. [TJ Red identifies the 
fits. The best fit values for 7 are provided for each contracted halo. The fits are generally 
fairly poor, although they improve at higher degrees of contraction and for higher a. 
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Fig. 3. — Fits to two different Einasto profiles at differing degrees of contraction, using 
double power laws (Eq. [Q. The profiles and colours are as found in Fig. [TJ Red identifies 
the fits. The best fit values for (3, 7 and S are provided for each contracted halo. The fits 
are considerably better than those found in Fig. |2j they tend to break down at higher a and 
at smaller degreees of contraction. 
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Fig. 4. — Scatter plot of the central power, 7 C , as a function of the shape parameter a, 
baryon fraction f^, baryon scale length and baryon to halo scale length ratio r^/r^ from 
left to right. The top panels refer to Blumenthal's AC prescription, while the bottom panels 
refer to Gnedin's AC prescription. The black points on all panels identify models for which 
0.1 < a < 0.2, which is consistent with cosmological simulations. The values for 7 C are largely 
constant with each of these quantities, and the aggregate values are given by 7 C = 0.75 ±0.17 
and 0.66 ± 0.08 for the Blumenthal and Gnedin prescriptions, respectively. Cases where 
r b > r h and r b > r outer are excluded from the plot. 
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Fig. 5. — The original halo rotation curve (black), contracted halo rotation curves (dotted 
red and blue) and the full (solid red and blue) AC rotation curves for the Blumenthal 
prescription (blue) and Gnedin prescription (red), at /b = 0.025, 0.05, 0.075, and 0.1 from 
top to bottom. Noise is added to these curves to create realistic mock data, but this is only 
shown for the highest rotation curve for clarity. The disc rotation curve for each model is 
identical and shown in green. 
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Fig. 6. — Fits to model rotation curve, adiabatically contracted for different values of /b, 
using plain Einasto profiles. Blue points with error bars identify some of the Blumenthal- 
contra cted mock data points a nd the blue curve identifies the best fit total rotation curve 
to the iBlumenthal et al.l (119861 ) prescription, while red points with error bars identify some 
of the Gnedin-contra cted mock data poin ts and the red curve identifies the best fit total 
rotation curve to the iGnedin et al.l (120041 ) prescription. For clarity, not all the mock data 
is shown. Green and orange identify the disc and halo breakdown of the rotation curves for 
the Blumenthal and Gnedin prescriptions, respectively. The disc rotation curve is shown as 
dashed lines, and the halo rotation curve is shown as dotted lines. The original disc (identical 
for all cases) is shown as the lowest black curve. The fits are very good, but overestimate 
the disc mass by 20 — 40%. 
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Fig. 7. — Fits to model rotation curve, adiabatically contracted for different values of /b, 
using NFW-type profiles in which (3 and 5 are held fixed. Colours and line types are as found 
in Fig. [61 The fits are similar to those in Fig. |6j 
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Fig. 8. — Fits to model rotation curve, adiabatically contracted for different values of /b, 
using double power laws in which all three indices (3, 7 and S are allowed to vary. Colours 
and line types are as found in Fig. El The fits are similar to those in Figs. [6] and [3 
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Fig. 9. — The rotation curve for NGC 6503, with the c ircula r vel ocity from PWC identified 
by blue circles and the ionized gas curve from lBottemal (119891 ) and lBegemanl (119871 ) identified 
by green triangles. 
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Fig. 10. — The r 2 p profile for NGC 6503, estimated from the circular velocity curve (top) and 
ionized gas curve (bottom). The circular velocity yields double humped profiles assuming 
no disc (black curve) and assuming a disc of 3.5 x 1O 9 M (green). The ionized gas curve 
provides noisy profiles with little useful data. 
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Fig. 11. — Fits to the rotation curve of NGC 6503 for the circular velocity fits. Black 
identifies the circular velocity curve obtained from PWC, plotted at radii where the ionized 
gas data was obtained. Red identifies the best fit rotation curve for fits that include a bulge, 
and orange identifies the rotation curves of the bulge (dotted), disc (dashed) and halo (solid) 
components. Blue identifies the best fit for fits that include no bulge, and green identifies the 
rotation curves of the disc (dashed) and halo (solid) components. The best fit x 2 f° r eacri 
scenario is indicated, for the bulge (top number) and bulgeless (bottom number) models. 
The labels B, G, and A correspond to the contraction prescriptions of Blumenthal, Gnedin, 
and Abadi, respectively. E corresponds to an uncontracted halo. 
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Fig. 12. — Fits to the r otation curve of NGC 6503 for t he ion ized gas fits. Black identifies 
the H/3/H I data from iBottemal (Il989l ) and iBegemanl (119871 ). Red identifies the best fit 
rotation curve for fits that include a bulge, and orange identifies the rotation curves of the 
bulge (dotted), disc (dashed) and halo (solid) components. Blue identifies the best fit for 
fits that include no bulge, and green identifies the rotation curves of the disc (dashed) and 
halo (solid) components. The best fit y 2 for each scenario is indicated, for the bulge (top 
number) and bulgeless (bottom number) models. The labels B, G, and A correspond to the 
contraction prescriptions of Blumenthal, Gnedin, and Abadi, respectively. E corresponds to 
an uncontr acted halo. 
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Fig. 13. — Stability parameters Q and X for each scenario. The colour pattern black-blue- 
red-green (going from top to bottom for each curve) identifies scenarios B, G, A and E, 
respectively. Each is plotted twice, once for the circular velocity fit (dot-dashed lines) and 
once for the ionized gas fit (solid lines). Only the profiles for models with a bulge are shown; 
profiles for the bulgeless models are similar. 
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Fig. 14. — Average X versus contraction R = Mf(rb)/Mi(rb) for each model that includes 
a bulge. Circles identify models from suite C and triangles identify models from suite I; 
the bulgeless models closely match these data points and are not plotted. Colours are as 
in Fig. [T3j The dotted line indicates the approximate cutoff for bar formation from PWC; 
models above the cutoff are less likely to display moderate bar formation. 
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Fig. 15. — PDFs for the Einasto shape parameter a. Blue identifies the PDFs for suite C, 
while green identifies the PDFs for suite I. Solid lines identify the fits that include a bulge, 
while dotted lines identify fits that are bulgeless. 
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Parameter 


Units 


Description 


Ph 


1O 9 M kpc- 3 


Halo density at r\ 




kpc 


Radius of r 2 p peak 


a 


Dimensionless 


Einasto shape parameter 


h 


Dimensionless 


Baryon fraction 




kpc 


Scale length of the final disc 


^outcr 


kpc 


Radius beyond which no baryons condense 


Calculated Quantities 


^disc 


1O 9 M 


Disc mass 


T^halo 


1O 9 M 


Halo mass inside r out er 


c 




Halo concentration 






Average X over two disc scale lengths 



Table 1: Table of the input parameters as well as calculated quantities that are output for 
each model. 
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Parameter 



Best-fit value 



Einasto fits, Blumenthal prescription 





/ = 


0.025 


}■ 


= 0.05 


/ = 


0.075 


/ 


= 0.1 


PO 


0.0022 


±0.0002 


0.0018 


±0.0002 


0.0018 


±0.0003 


0.0020 


±0.0004 


rh 


17.6 


±0.6 


19.0 


±1.3 


18.4 


±1.2 


17.3 


±1.6 


a 


0.130 


±0.003 


0.123 


±0.004 


0.125 


±0.004 


0.128 


±0.006 


^disc 


40.3 


±1.7 


40.7 


±1.4 


39.8 


±1.1 


38.8 


±1.1 


Einasto fits, Gnedin prescription 




/ = 


0.025 


/ = 


= 0.05 


/ = 


0.075 


/ 


= 0.1 


PO 


0.0023 


±0.0002 


0.0018 


±0.0002 


0.0018 


±0.0003 


0.0019 


±0.0003 




17.1 


±0.6 


19.2 


±1.1 


19.0 


±1.6 


17.8 


±1.5 


a 


0.125 


±0.003 


0.118 


±0.004 


0.119 


±0.005 


0.123 


±0.005 


m disc 


34.3 


±1.6 


37.8 


±1.3 


37.9 


±1.3 


37.4 


±1.0 


NFW fits, Blumenthal prescription 




/ = 


0.025 


/ = 


= 0.05 


/ = 


0.075 


/ 


= 0.1 


Ps 


0.0017 


±0.0002 


0.0017 


±0.0002 


0.0018 


±0.0003 


0.0024 


±0.0005 


r s 


34.6 


±1.4 


33.6 


±2.0 


31.7 


±2.3 


27.1 


±2.5 


a 


3 




3 




3 




3 






1 




1 




1 




1 




7 


1.35 


±0.02 


1.36 


±0.01 


1.35 


±0.01 


1.33 


±0.01 


m disc 


49.5 


±1.6 


47.6 


±1.3 


46.3 


±1.1 


44.2 


±1.1 


NFW fits, Gnedin prescription 




/ = 


0.025 


/ = 


= 0.05 


/ = 


0.075 


/ 


= 0.1 


Ps 


0.0015 


±0.0001 


0.0015 


±0.0002 


0.0016 


±0.0002 


0.0032 


±0.0005 


r s 


36.4 


±1.4 


35.8 


±2.0 


34.3 


±2.3 


24.0 


±1.6 


a 


3 




3 




3 




3 







1 




1 




1 




1 




7 


1.38 


±0.02 


1.38 


±0.01 


1.37 


±0.01 


1.33 


±0.01 


^disc 


44.2 


±1.6 


44.9 


±1.2 


44.9 


±1.0 


41.2 


±1.0 


Double power fits, Blumenthal prescription 




/ = 


0.025 


/ = 


= 0.05 


/ = 


0.075 


/ 


= 0.1 


p s 


0.0042 


±0.0004 


0.0037 


±0.0005 


0.0062 


±0.0003 


0.0068 


±0.0003 




28.7 


±3.8 


23.5 


±3.2 


40.4 


±2.3 


30.9 


±2.2 


a 


3.01 


±0.12 


2.76 


±0.13 


3.30 


±0.09 


3.15 


±0.13 




0.76 


±0.09 


0.92 


±0.16 


0.52 


±0.02 


0.56 


±0.03 


7 


1.24 


±0.04 


1.31 


±0.04 


1.11 


±0.02 


1.13 


±0.02 


m disc 


44.9 


±1.7 


46.0 


±1.5 


42.4 


±1.1 


41.5 


±1.3 






Double power 


fits, Gnedin prescription 








/ = 


0.025 


/ = 


= 0.05 


/ = 


0.075 


/ 


= 0.1 




0.0044 


±0.0004 


0.0044 


±0.0002 


0.0069 


±0.0005 


0.0068 


±0.0006 


r s 


35.0 


±1.4 


39.3 


±3.1 


37.3 


±1.7 


32.3 


±2.2 


a 


3.17 


±0.04 


3.18 


±0.10 


3.22 


±0.11 


3.17 


±0.12 




0.64 


±0.02 


0.57 


±0.03 


0.52 


±0.02 


0.56 


±0.03 


7 


1.21 


±0.02 


1.19 


±0.02 


1.13 


±0.01 


1.14 


±0.01 


m disc 


38.1 


±1.6 


40.5 


±1.3 


40.1 


±1.5 


39.6 


±1.4 



Table 2: Posterior mean values for the input parameters with 1— a error bars for the suite of 
test rotation curve fits. 
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Inner rotation fits, Blumenthal prescription 



/ = 0.025 / = 0.05 / = 0.075 7 = 0.1 





0.63 


±0.18 


0.91 


±0.28 


0.96 


±0.18 


0.75 


±0.21 


7 


0.77 


±0.09 


0.92 


±0.09 


0.93 


±0.07 


0.80 


±0.08 


5 


3.80 


±1.07 


6.92 


±1.27 


6.85 


±1.27 


4.03 


±0.81 






Inner rotation fits, Gnedin prescript 


ion 






/ = 


0.025 


/ = 


0.05 


/ = 


0.075 


/ 


= 0.1 


P 


0.50 


±0.08 


1.02 


±0.57 


0.79 


±0.23 


0.90 


±0.19 


7 


0.67 


±0.05 


0.87 


±0.13 


0.82 


±0.09 


0.88 


±0.07 


5 


5.11 


±1.07 


4.37 


±1.28 


4.03 


±0.92 


4.59 


±0.79 



Table 3: Fitted double power shape parameters to the inner (<5 kpc) rotation curves of the 
AC halos in §2.2, with 1 — a error bars. 



Shape parameters, Blumenthal prescription 





f = 0.025 


/ = 0.05 


/ = 0.075 


/ = 0.1 




0.59 


0.53 


0.53 


0.56 


7 


0.76 


0.76 


0.76 


0.77 


S 


3.23 


4.95 


4.84 


5.04 



Shape parameters, Gnedin prescription 
/ = 0.025 / = 0.05 / = 0.075 / = 0.1 





0.50 


0.55 


0.54 


0.54 


7 


0.60 


0.63 


0.64 


0.64 


5 


3.40 


3.13 


3.08 


3.03 



Table 4: Double power shape parameters for the inner profiles of the AC halos in §2.2, 
obtained directly from the r 2 p profiles. 
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Parameter 



Best-fit value 



Suite I 







B 




G 




A 




E 


Ph 


0.0024 


±0.0003 


0.0025 


±0.0004 


0.0024 


±0.0004 


0.0029 


±0.0004 




11.8 


±0.9 


11.6 


±1.1 


11.6 


±0.9 


10.1 


±0.6 


a 


0.853 


±0.073 


0.885 


±0.061 


0.664 


±0.186 


0.729 


±0.142 


rb 


1.34 


±0.04 


1.40 


±0.04 


1.22 


±0.04 


1.3 


±0.05 


fb 


0.058 


±0.005 


0.076 


±0.005 


0.064 


±0.016 






Pouter 


71.6 


±22.3 


72.7 


±20.8 


70.4 


±19.2 






™halo 


100 


±10 


95 


±8 


114 


±21 


90 


±18 


m-disk 


5.8 


±0.4 


7.2 


±0.5 


7.0 


±0.9 


9.2 


±0.8 


C 


7.91 


±0.44 


7.95 


±0.51 


8.41 


±0.60 


8.56 


±0.53 


(X) 


1.50 


±0.06 


1.28 


±0.05 


1.06 


±0.09 


0.92 


±0.04 


x 2 


2.78 




2.69 




2.45 




2.38 




Suite I, no bulge 






B 




G 




A 




E 


Ph 


0.0029 


±0.0003 


0.0030 


±0.0003 


0.0031 


±0.0003 


0.0036 


±0.00038 


rh 


10.8 


±0.6 


10.6 


±0.5 


10.4 


±0.6 


8.9 


±0.5 


a 


0.850 


±0.086 


0.901 


±0.051 


0.769 


±0.105 


0.541 


±0.139 


rb 


1.24 


±0.03 


1.29 


±0.04 


1.11 


±0.04 


1.17 


±0.05 


fb 


0.057 


±0.005 


0.077 


±0.005 


0.078 


±0.011 






Pouter 


72.6 


±20.2 


68.9 


±21.6 


47.0 


±21.6 






"Ihalo 


94 


±8 


87 


±5 


88 


±10 


101 


±17 


m-disk 


5.3 


±0.3 


6.7 


±0.6 


6.8 


±0.6 


7.4 


±0.8 


C 


8.41 


±0.33 


8.39 


±0.32 


8.57 


±0.57 


10.0 


±0.9 


(X) 


1.42 


±0.06 


1.21 


±0.05 


0.96 


±0.04 


0.95 


±0.06 


x 2 


1.87 




1.82 




1.75 




1.75 




Suite C 






B 




G 




A 




E 


Ph 


0.0028 


±0.0001 


0.0028 


±0.0001 


0.0028 


±0.0001 


0.0036 


±0.0002 


rb 


9.9 


±0.2 


9.7 


±0.2 


9.7 


±0.1 


8.5 


±0.2 


a 


0.120 


±0.005 


0.133 


±0.005 


0.107 


±0.003 


0.101 


±0.003 


rb 


1.38 


±0.02 


1.41 


±0.02 


1.33 


±0.02 


1.35 


±0.02 


fb 


0.026 


±0.003 


0.029 


±0.002 


0.026 


±0.002 






Pouter 


38.8 


±5.7 


41.1 


±5.1 


45.5 


±4.3 






W-halo 


117 


±14 


121 


±12 


134 


±10 


214 


±11 


m-disk 


3.0 


±0.1 


3.4 


±0.1 


3.5 


±0.1 


4.3 


±0.1 


C 


9.89 


±0.41 


10.20 


±0.28 


10.60 


±0.32 


10.20 


±0.91 


(X) 


3.68 


±0.10 


3.27 


±0.09 


2.95 


±0.07 


2.44 


±0.07 


x 2 


2.41 




2.49 




2.24 




2.06 




Suite C, no bulge 






B 




G 




A 




E 


Ph 


0.0027 


±0.0002 


0.0024 


±0.0002 


0.0028 


±0.0001 


0.0027 


±0.0002 


rb 


9.9 


±0.3 


10.6 


±0.5 


9.8 


±0.2 


9.7 


±0.3 


a 


0.093 


±0.003 


0.098 


±0.004 


0.087 


±0.003 


0.078 


±0.002 


rb 


1.42 


±0.02 


1.46 


±0.02 


1.36 


±0.02 


1.39 


±0.02 


fb 


0.024 


±0.003 


0.032 


±0.005 


0.027 


±0.002 






Pouter 


37.4 


±4.2 


33.6 


±5.3 


40.2 


±3.6 






"*halo 


115 


±11 


104 


±14 


123 


±9 


242 


±19 


m-disk 


2.8 


±0.1 


3.2 


±0.1 


3.3 


±0.1 


4.4 


±0.1 


C 


9.83 


±0.54 


8.93 


±0.73 


10.20 


±0.46 


9.39 


±1.25 


(X) 


4.04 


±0.10 


3.57 


±0.09 


3.20 


±0.08 


2.51 


±0.06 


x 2 


2.56 




2.72 




2.36 




2.25 





Table 5: Table of the best fit parameter values for NGC 6503 rotation curve fits with 1— a 
error bars. 
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